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A nonreflecting wavepacket is constructed by the superposition of reflectionless eigenstates of 
sech2 potential. Free propagation and propagation in the presence of the above potential of such 
a wavepacket is considered using the concept of arrival time. Comparison is made with the free 
evolving Gaussian wavepacket. Mean arrival time at a detector behind the well is given as a function 
of mass for separate cases. A selection of Bohmian trajectories in the interacting case are computed 
and compared to the trajectories of a free particle guided by a Gaussian packet. 
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I. INTRODUCTION 

Scattering of incident waves by inhomogeneities is a common physical feature in optics and quantum mechanics. 
The trigonometric Poschl-Teller (PT) potential was suggested by Poschl and Teller [ij to describe diatomic molecular 
vibration. PT-type potentials have been noticed by authors from different areas of physics ranging from the molecular 
physics to the nuclear one. Using different approximation methods, analytical solutions of the Schrodinger equation 
and energy eigenvalues for diatomic molecules have been found for the trigonometric PT molecular potential 0. 
Approximate bound-state solutions of the Klein-Gordon Q and Dirac equations with such potentials have been 
investigated. The Sech-squared potential (or Modihed PT potential) has attracted too too much attention, for instance 
it appears in the context of neutron-matter Q and spherical quantum dot [^. 

There are obstacles, called reflectionless potentials, which are completely transparent for the incident wave. This 
type of potentials have been extensively studied analytically and numerically in the literature. Using the method 
of summation over eigenstates, the exact propagator for a general reflectionless potential has been computed Q- 
Scattering of a Gaussian wavepacket by a reflectionless sech2 potential has been numerically studied and possible 
applications of such potentials has been argued 0 . It has been shown that in comparison to the case of free evolution, 
the wavepacket accelerates in the well and is narrower. Authors arrived at this result by plotting the probability 
density versus position in a time that the whole wavepacket has left the well. In this paper we employ the concept 
of arrival time Q for quantitative study of this problem. In the case of reflectionless sech2 potential non-reflecting 
wavepackets have been constructed by superposing energy eigenstates and propagation of them has been studied in 
the presence of the above potential [^. Similar results have been reproduced by the method of supersymmetric 
quantum mechanics llll . Traversal time through the reflectionless sech2 potential has been studied using the method 
of a quantum clock |12l |. Only, recently, a true reflectionless potential was realized experimentally in photonic lattices 
using the concept of arrays of evanescently couple waveguides [l^ . 

It seems to be interesting to consider the problem of a wavepacket scattering by the reflectionless Sech2 potential 
using a trajectory-based theory, since by attributing a well-dehned trajectory to a particle, its arrival time at a 
detector location is a well-defined and predictable quantity. It has been shown within the Bohmian mechanics [l^ 
that the arrival time distribution is given by the probability current density [Tsj l. But, there is no consensus about 
this quantity within the standard formulation of the quantum mechanics [l6j |. 

To this end, we at first follow Lekner 0 to construct a non-reflecting wavepacket. Then, propagation of such a 
constructed wavepacket is considered and arrival time at a detector location behind the well is computed. Gomparison 
with a freely evolving Gaussian wavepacket is made. Finally, a selection of Bohmian paths is plotted. 
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A. Constructing a nonreflecting wavepacket 

^2i2 . . 

Exact solution of the Schrodinger equation for any positive energy is given by 
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in the presence of the sech2 potential 
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with V = \. One easily sees that in the limit x —>• —oo, the above solution reduces to A(k)e^^^. This shows the 
nonreflectivity of the potential; with an incident wave from the left there is no reflected wave [l7| . 

By superposing the reflectionless eigenstates o, one obtains nonreflecting wavepacket 
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where sub-indices ”i” and ”nr” stand respectively for ’’interacting” and ’’nonreflecting”. Noting the factor (ik + a) ^ 
in the relation of il^kix), a simple choice for A{k) is 0 
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where, 
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is the Fourier transform of the initial Gaussian wavepacket 
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where subindex ” G” stands for ” Gaussian”. Here, /cq is the kick momentum; and Xc and cto are respectively the center 
and rms width of the packet. Time-evolved Gaussian wavepacket in free space is given by. 
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where St = cto(1 +ifit/2maQ) and u = fik^/m is the group velocity of the packet. In this equation subindex ”f’ stands 
for ” free”. 

By the choice (|1]) for H(fc), reflectionless wavepacket reads 
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where iV, which has the length dimension, shows the normalization constant. 

Free evolution of the initial nonreflecting wavepacket = dij nr(a;, 0) is given by the integral relation 
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where. 
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represents the propagator of free particle. From eqs. m and one obtains 
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where, 


g{x,t) = / dx' tanh(aa:') exp <— 
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Just for completeness we mention that the evolution of the Gaussian wavepacket ([6|) in the presence of the nonre¬ 
flecting potential ([2) with i/ = 1 which has been numerically studied in is given by 
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represents the propagator of the nonreflecting potential for :/ = 1 


II. BASIC EQUATIONS 

Using eq. ([8|) probability density and probability current density becomes 


Pi,nr{x,t) = |fV|" 
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where at = |st| = cro^/l-j- and 
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According to the Ehrenfest theorem, expectation value of momentum remains unchanged in the free space. Thus, 
in the free space one has 


{x){t) = (a;)(0) + MMt. 

m 

For the freely evolving Gaussian wavepacket this equation reads, 
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In the Bohmian view of the world, the state of a system is described by its wave function together with its actual 
configuration defined by the actual position of its particle. In the minimal approach to the theory [T^ . particle 
equation of motion, known as guidance equation, is 
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Using this equation and eqs. m and (US, one obtains 
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for the initial Bohmian velocity of a particle which is guided by the nonreflecting wave packet It is clear from eq. 
(El) that initial velocity depends on the position in the nonreflecting wavepacket case, while it is constant value u for 
the Gaussian packet. By integrating guidance equation (1161) over time one obtains 


a:f.G(i) = Xc + ut +—{xq - Xc) , (18) 

for the trajectories of a free particle which is guided by the Gaussian packet ©■ Here, xq is the initial position of the 
particle. 

Within the Bohmian mechanics, the arrival time distribution of particles at detector location Xd is given by flB | 


and thus the mean arrival time reads 
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III. SCALING THE UNIT OF LENGTH 


In order to do numerical calculations, we apply a scaling transformation which eliminates physical constant fi 
and mass m in the Schrddinger equation for free particle and for a particle interacting with sech2 potential ([2]). 
Wavefunction 4'(a;,t) which describes a particle of mass m in the potential ([2]) with 1 ^ = 1, satisfies the Schrddinger 
equation 
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By space transformation x —>■ x = y/Wx, which is a scaling of unit length [20|, and defining new function 
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Schrddinger equation (EU reads 
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where a = y ^a. Thus, transformed momentum operator is p = —idfdx. 
transformation corresponds to put = to = 1 in the equations of section im 
Gontinuity equation becomes 
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says that the above scale 
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and thus, guidance equation in this new coordinate reads 
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In summary we have the following transformations, 
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IV. NUMERICAL CALCULATIONS 

For numerical calculations, transformed parameters of the potential and wavepacket are taken as a = 1, do = 1; 
Xc = —10, kg = 1. In figure [U we have plotted transformed probability density p(x,t) versus transformed space 
coordinate x at six different instants. A gap at the center of the potential is seen in the plot of interacting nonreflecting 
wavepacket. With respect to the relation Xd = yjm/h Xd, for a given width 6x of p, corresponding width in x is 
proportional to the inverse of square root of mass. Thus, from the forth equation of eq. it follows that by increasing 
the mass, the width of probability density decreases . With the above values for the parameters of the wavepacket 
one obtains by numerical integration N ~ 2.25 for the normalization constant in eq. (p)nr(O) = (^nr|p|^nr)(0) 
1.2222 and (a:)jjj.(0) = ('I'nr|a;|4'nr)(0) — —10.4444 for the initial nonreflecting wavepacket. Because of the relation 
(p)jjj.(0) > (p)g(O), expectation value of space coordinate of the nonreflecting wavepacket, (x)^^, in the free-space 
moves faster than the center of the Gaussian wavepacket. In addition, in the interacting case because potential is 
attractive, wavepacket moves faster when the wavepacket approaches the well. After, passing the potential it moves 
with constat velocity. See Fig. [21 It must be noted that potential is proportional to the inverse of mass. Thus, it is 
more attractive for low masses. So, wavepackets describing particles with low values of mass move faster in comparison 
to the large ones. As one sees in plot b) of figure H] uncertainty in position, Ax = ^(tknrlx^ldlnr) — (4'nr|a;|'I'nr)^, is 
larger for the interacting nonreflecting wavepacket compared to the freely evolving Gaussian wavepacket. 

In figure [3] we have plotted arrival time distribution (1191) at the detector location Xd = 2d. Relation Xd = \/m/fi Xd 
means fixing Xd and changing Xd corresponds to changing mass. From plots (a) and (b) of this figure it is found that 
the width-in-time of the arrival time distribution increases with Xd (mass) in contrast to the width-in-space of the 
detection probability density p where width-in-space in a given time decreases with mass. Plot (c) of Fig. [3] shows 
the mean arrival time (1201) . In agreement with the previous findings, f is the lowest for the interacting nonreflecting 
wavepacket while is the largest for the free Gaussian. For all types of packets mean arrival time increases with Xd 
(mass). 

Figure m represents a selection of Bohmian paths. In plot (a), trajectories (TTSl) for the free Gaussian have been 
plotted, while in plot (c) trajectories are for a particle which is guided by the wavepacket dt, nj.(a:, t). In agreement 
with plot (c) of Fig. O trajectories spread out more in the case of interacting nonreflecting packet. In plot (c) we 
have depicted {x){t) and the Bohm path with initial position (a;)jjj,(0). It is seen that the Bohm path overtakes the 
expectation value of space coordinate. Plot (d) represents initial Bohmian velocity (TTfl) versus space coordinate x. 
As one sees, velocity is higher for the interacting nonreflecting wavepacket compared to the free Gaussian. 
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FIG. 1: (Color online) Transformed probability density p(x, t) versus transformed space coordinate x at six different instants 
t = 0, t = 4 s, t = 8 s, t = 12 s, t = 16 s and t = 20 s. In each plot black curve shows interacting nonreflecting wavepacket while 
the red one shows freely evolving nonreflecting packet and the green one stands for the propagation of free Gaussian packet. 



FIG. 2: (Golor online) Expectation value of space coordinate and momentum (a and b); and uncertainty in position (c) with 
the same color scheme as in Fig. [T] 
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FIG. 3: (Color online) Arrival time distribution Il{xd, i) versus time (left column) and mean arrival time T{xd) versus detector 
position Xd with the same color scheme as in Fig. [T] 

V. SUMMARY AND CONCLUSION 

In this work we constructed a nonreflecting wavepacket by superposing of eigenstates of the reflectionless sech2 
potential. Then, evolution of this wavepacket in the presence of the potential and in free space studied and compared 
to the free propagation of Gaussian wavepacket. It was found that nonreflecting wavepacket moves faster and its 
width, uncertainty in position, is larger. Mean arrival time at a detector behind the well was computed as a function 
of mass. It was seen mean arrival time increases with mass in agreement with one’s intuition and it is the largest 
for the free Gaussian compared to the nonreflecting wavepacket in free space and in the interacting case. Finally, a 
selection of Bohmian trajectories were computed. 
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